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ABSTRACT 

We examine the constraints imposed by helioseismic data on the solar heavy el- 
ement abundances. In prior work we argued that the measured depth of the surface 
convection zone Rcz and the surface helium abundance Y^urf were good metallicity in- 
dicators which placed separable constraints on light metals (CNONe) and the heavier 
species with good relative meteoritic abundances. The resulting interiors-based abun- 
dance scale was higher than some published studies based on 3D model atmospheres at 
a highly significant level. In this paper we explore the usage of the solar sound speed in 
the radiative interior as an additional diagnostic, and find that it is sensitive to changes 
in the Ne/0 ratio even for models constructed to have the same Rcz and Ysurf- Three 
distinct helioseismic tests (opacity in the radiative core, ionization in the convection 
zone, and the core mean molecular weight) yield consistent results. Our preferred O, 
Ne and Fe abundances are 8.86 ± 0.04, 8. 15 ± 0. 17 and 7.50 ± 0.05 respectively. They 
are consistent with the midrange of recently published 3D atmospheric abundances 
measurements. The values for O, Ne and Fe which combine interiors and atmospheric 
inferences are 8.83 ± 0.04, 8.08 ± 0.09 and 7.49 ± 0.04 respectively. 

Subject headings: atomic data - opacity - Sun: abundances - Sun: evolution - Sun: interior 
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1. Introduction 



The Sun is a fundamental benchmark in astronomical abundance studies. Meteorites can 
be used to reliably infer relative abundances for numerous heavy elements and we have detailed 
information about limb darkening, precise fundamental parameters, and a wide variety of 
diagnostic features. A decade ago the problem of solar abundances seemed settled. A series 



of papers (I Anders & Grevesse 



1989 



Grevesse & Noels 



1993 



Grevesse & Sauval 



1998 



, GS98 



hereafter) had estab lished a consistent set of surface abundances, and interiors models (e.g. 



Bahcall et al. 



200 1|) were in good agreement with helioseismic data ranging from sound speed as 



a function of depth to surface helium and convection zone depth measurements. However, the 
atmospheres models used in this generation of study were one-dimensional and had important 
simplifying assumptions about the physics of the atmospheres. With the advent of a new 
generation of numerical simulations it became possible to relax these severe assumptions. 
For many heavier elements such as iron the revisions in abundance were small and within 
the errors. However, the lighter elements (CNO) were a different story. These species have 
problematic photospheric abundance diagnostics under solar conditions and meteorites cannot 
be used to infer their primordial relative abundGances. The 3D model atmospheres, along with 
other revisions such as enhanced non-LTE effects and the identification of blending features in 



commonly used abundance indicators, led to much lower CNO abundances (|Asplund et al 



2005 



AGS05 hereafter). The revised abundanc es placed solar interiors models in stark disagreeme nt 



with helioseismic data (see, for examp 



Bahcall et al. 



2005al : 



Guzik et al. 



2005 



Antia & Basu 



Montalban et al. 



2005 



2004; 



Bahcall & Pinsonneault 



Turck-Chieze et al 



solar abundance pr oblem spurred a re-examination of the opacities (|Seaton & Badnell 



Badnell et al 



Opacity Project (jSeaton et al. 



2004 



2004. Th is 



2OO4I : 



20051). which improved opacity data as compa red to the original source from the 



19941 : 



review by Zeippen 



1995h but did not resolve the problem. 



The next generation of helioseismic studies took a different approach, inverting the problem 



and u sing distinct characteristics to infer a higher seismic abundance more in li ne with the GS98 



scale ( Antia & Basu 



2007 



Chaplin et al. 



2007 



Delahave & Pinsonneault 



2006L DP06 hereafter). 



DP06 identified the depth of the surface convection zone and the surface helium abundance as 
sensitive abundance diagnostics. This data could be used to predict absolute light and heavy 
metal abundances when combined with relative meteo ritic abundances for heavy metals and 
relative photospheric/coronal abundances for CNONe. .Antia & Basul (|2007h used the depression 
of the adiabatic gradient induced by metal ionization, mostly from oxygen, to infer the absolute 



metal 



icty of the solar convection zone and obtained a comparably high abundance. 



Chaplin et al 



(l2007h derived a similar result from modelling the metallicity dependance of the mean molecular 



weight in the deep solar core. 

There was a subsequent burst of revisions in the atmospheric abundances. The initial 
generation of models were criticised for problems in fitting the observed solar limb-darkening 



(lAyres et al 



20061 : 



Koesterke et al 



20081). Substantial revisions of the original model atmospheres 



led to a signficant improvement in the fit to solar data and remov ed the mild tension b etween the 



abundances of heavier elements (Fe, Si) and seismic inferences (|Asplund et al. 



2009L hereafter 



AGSS09). However, the estimates for CNO abundances remained quite low. Interestingly, 
a different group of investigators with an independant 3D atmospheres code had compara ble 



agree ment with solar data but found higher photo spheric abundances (|Caffau et al. 



2008L 



2010|) . The solar abundance problem therefore takes on a different cast; the distinction is 
not between new 3D and old ID atmospheres, but rather on differing choices of abundance 
indicators and different ways of analyzing them. A similar point was made by us in a prior paper 



(IPinsonneault & Delahave 



2009|) . where a reanalysis of the AGS05 work found an atmospheric 



abundance scale intermediate b etween seismology and the prior 3D values, comparable to the one 



obtained by 



Caffau et al. 



( 2010 ) 



However, neon has always been an important, and underconstrained, ingredient in the puzzle. 



Ne is a significant opacity source in stellar interiors, so the helioseismic pr oblems induced by 



low oxygen could be compensated for if the solar Ne was sufficiently high (lAntia & Basu 



2005 



iahcalletal 



( Turck-Chieze et al. 



2005b). The re are also chemical evolution arguments in favor of such a model 



2004) . Ne is not preserved in meteorites in significant amounts and lacks 



photospheric diagnostics because of its atomic structure. As a result, indirect arguments from 



coronal abundance ratios must be usee 



been argued for (|Drake & Testa 



2005 



, and significant variation s from the standard mixture have 



Drake & Ercolano 



2005 



2007). However, the bulk of evidence 



Young 



2005); see AGSS09 for a 



in the Sun appears to favor a lower ratio (|Schmelz et al 
detailed discussion. 

In DP06 we used only the scalar constraints (Rcz and Fvur/) as diagnostic features to infer 
abundances. However, Ne retains electrons to higher temperatures than O does, and as a result 
different Ne+0 mixtures with the same opacity at the base of the surface convection zone will 
have different opacity in deeper layers. This leads to a distinct signature in the sound speed profile. 
In this letter we update the mixture constraints based on the same method, but with improved 
stellar interiors model physics. By adding the sound speed profile as an additional constraint we 
can set interesting bounds on the solar Ne/0 ratio, and these are reinforced by considering the 
seismic metallicity bounds based on ionization and mean molecular weight. In section 2 we recall 
the method used in DP06, to determine the appropriate composition in order to reproduce the 
two observable scalars Rcz and Y^urf and the extension of this method including the additional 
constraint Csound to determine the ratio Ne/0. In Section 3 we discuss our results. 

2. Method 



In DP06 we presented our method to infer the solar heavy and light metal abundances 
required for theoretical models that reproduce the observed solar Ysurf and Rcz- There are three 
essential steps involved in determining a seismic solar composition and its associated error. 
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First, the observed solar properties and their errors must be obtained. Second, errors in the 
input solar model physics introduce uncertainties in the theoretically predicted K,„r/ and Rcz- 
These uncertainties can be correlated; for example, increasing the degree of gravitational settling 
both deepens the model Rcz and decreases the K,„,-/. Both the difference between theory and 
observation and the associated error are now defined for a given solar composition. 

Changes in the abundance of the heavier metals primarily affect the Y^urf through their impact 
on the central opacity and initial helium, while changes in the CNO abundances primarily affect 
Rcz ■ As a result, we can solve for the heavy (relative abundances inferred from meteorites) and 
light (relative abundances inferred from photospheric and coronal data) abundances consistent 
with the solar data as a function of the assumed Ne/0. We derived the seismic solar abundances 
in term of A(Fe) and A{0) keeping the ratios C /O, N/0, Ne/0 and Xi/Fe (where Xj stands for all 
metals but CNONe) constant. A{X) is the abundance of element X on the logarithmic scale with 
A{H) = 12. We obtained A(C>)ag505 = 8.87 dex±0.041 and A(Fe)AGS05 =7.51 dex±0M5 for the 
calibrated model using the relative CNONe abundances from the AGS05 compilation. Similar 
results were obtained using relative light element abundances from GS98. 

For the present work w e first updated our solar model inputs, using an equation of state table 



([Rogers & Nayfonovl2002|) with revised online data (1http://adg.llnl.gov/Research/OPAL/opal.html 



hereafter EOS2006), l ow temperature opacity tables from 



nucle ar reaction rates (|Junghans et al 



2003 



Brown et al, 



Ferguson et al. 



2007 



Confortola et al 



(l2005h and updated 



2007 



Martaetal 



20081). We then solve for the uncertainties in sound speed profile for models with fixed Rcz 
and Ysurf- This defines a range of sound speed deviations which cannot be distinguished from 



theoretical errors, and which therefore set a bound on our ability to infer the Ne/0 ratio. 
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2.1. Uncertainty envelope for the sound speed profile 



The uncertainties associated with the sound speed profile icsound) include the observational 
errors in c^^J^^, uncertainties in the predicted Csoimd due to the uncertainties in Rcz and Y^urf and 
finally the theoretical errors from uncertainties in the model input physics (opacities, equation of 
state, nuclear reaction cross sections, rate of gravitational settling) and assumed solar properties 
(age, radius, and luminosity). In the top panel of Figure 1 we illustrate individual theoretical 
uncertainties with lines and the total effect by the gray shade. On the bottom panel the measured 
Csound uncertainties are in light gray, the ones induced by observational uncertainties in Rcz and 
Ysurf are in medium light gray and the theoretical error envelope is in medium dark gray. The total 
error envelope is represented in dark gray in the bottom panel. Our detailed procedure is described 
below. 

In order to estimate the effect of the different physics ingredients onto the Csound profile, we 
changed one at a time and adjusted the composition in the same manner as in DP06 to recover the 
seismic Rcz and Ysm-f- The differences in Csomd between these calibrated models are then measures 
of the impact of theoretical uncertainties on the Csound profile. We adopt the same error estimates 
and central values as described in Table I of DP06 with the following exceptions. The errors 
and central values of the nuclear reaction cross sections have been taken from the sources cited 
above and for the equation of state the errors are derived from the differences in the prediction 
between the EOS2001 and EOS2006 tables. We used the differences in the calibrated model 
sound speed profiles to infer the envelope of uncertainties. From the top panel of Figure 1, we 
see the uncertainties in the diffusion coefficients for gravitational settling (dot-dashed line) are the 
most significant error source, followed by the uncertainties in the opacities (solid line). 

The observat ional uncertainties in the solar csound are deduced from the range of inverted 



Basuetal 



(|2000|) (see Figure 1, bottom panel, light gray envelope). Errors 



solar structures in 

in Rcz and Y^urf also impact the csound profile. The fractional depth of the convection zone is 
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Rcz = 0.7133 ±0.0005 Rq (iBasu & Antia 



2004). The theoretical location of the convection zone 



depth can change by up to 0.001 with different opacity interpolation methods (IBahcall et al. 



20041); we therefore treat the latter as the uncertainty in Rcz- The helium abundance in the surface 
convection zone {Ysurf = 0.2483 ± 0.0046) is me asured frorn the strength of the depression in the 



1998 



Richard et al. 



1998 



Di Mauro et al. 



adiab atic gradient in the helium ionization zone (|Basu 
2OO2L see DP06 for discussion). The observational errors in Rcz and Y^urf permit a range of 
derived compositions which will then generates different predicted Csound profiles. The values 
Rcz'" and Y^'^fjf define a reference mixture. The maximum Csound profile change occurs for 
^'cz =Rcz^" + (yRHeiw and y/„^y = Y^fr" + aYHeiw. This is our one a errors envelope (the medium light 
gray shade in Figure 1, bottom panel). We added the three sources defined above in quadrature to 
obtain the total envelope for uncertainties in the sound speed profile difference (dark grey shade 
of Figure 1 bottom panel). The theoretical errors predominate. 



3. Constraining Ne/0 



We are now in a position to extract additional information from the sound speed profile 
above and beyond those constraints imposed by the measurement of Rcz and K,,,,-/- Carbon 
and nitrogen are much less abundant than oxygen and have similar ionization structures. As a 
result, our method does not place interesting constraints on the C/0 and N /O ratios. However 
changes in Ne/0 have a significant impact on the csound profile, as seen in the top panel of 
Figure 2 where we compare different calibrated models to one reference model which has 
(0,Ne,Fe,Ne/0) = (8.80,8.25,7.49,0.28). In the bottom panel we present the percentage 
difference between the solar Csound and the different calibrated models with A(0) ranging from 
8.90 (top) to 8.70 (bottom) in 0.05 dex increaments. The dashed line correspond to our prefered 
model which minimizes the rms in the percentage difference in csound- Ne remains opaque to 
higher temperatures than oxygen, so a combination of higher Ne and lower O which yields the 



same opacity at the base of the convection zone will be more opaque deeper in the star. This 
difference is at a maximum relative to the theoretical uncertainties at 0.56 solar radii. Changes 
of A(A'^e) of 0.17 dex at this radius are comparable to the 1 a total error. We treat 0.17 dex as 
our uncertainty on the Ne/0 ratio from the c sound profile independent of the scalar or coronal 
abundance constraints. 

In Figure 3 we present the predicted Rcz and Ysurf (top panel) and the percentage difference 
in Csoimd between the Sun and severals models using various compositions. As previously noted, 
low O scales are strongly disfavored in both the scalar constraints (top panel) and the Csound profile 
(bottom panel). All of the models are discrepant at a statistically significant level in the CsouM 
immediately below the surface convection zo ne. We attribiite this to the treatment of the mixing 



which is required to explain the low solar Li (IPinsonneauh 



19971) . As in DP06, we reduced the 



settling coefficient to account for the eff ect of rotational mixing from full evolutionary models 



(IRichard et al 



19961 : 



Bahcall et al. 



20011). However rotational mixing requires a more complete 



model so we don't use this feature as a diagnostic and defer such a topic to a subsequend paper. 
Also illustrated in the bottom panel of Figure 3 (dashed line) is a model with low O (value from 
AGSS09) and high Ne (to satisfy the helioseismic contraints Rcz and Ysurf) which is excluded at 
more than two a. 



4. Conclusion 

Helioseismology has proven to be a powerful means of inferring the solar composition. 
The depth of the surface convection zone and the initial helium abundance are very sensitive 
to the opacity of light and heavy metal respectively, and the errors in these model properties 
can be reliably quantified. Reproducing these seismic scalar features require a relatively high 
solar metallicity. If we combine these scalar constraints with reasonable priors about the 
relative abundances of CNONe (from the photosphere) and heavy elements (from meteorites) 
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we can therefore infer a seismic solar mixture based primarily on stellar interiors physics. 
Recent revisions to the input physics of the solar model - improvements in the nuclear reaction 
cross-sections, low temperature opacities, and equation of state - yield a reference A(0) = 8.86 
and A(Fe) = 7.50 which is the same as DP06. 

However, the linkage between the thermal structure and the composition is indirect, which 
can lead to ambiguities in the interpretation of the strong constraints imposed by the measured 
solar thermal structure. Fortunately, there is additional information encoded in the solar sound 
speed profile and there are independent helioseismic diagnostics which yield very similar results. 

It is common to use differences in sound speed between models and the Sun as a measure 
of agreement. However, many of the striking visual features in such difference plots are in 
fact manifestations of errors in the independently determined scalar features, especially the 
convection zone depth. We therefore argue that extracting information from the sound speed 
profile first requires calibration to a common set of observables, which can be achieved with 
modest adjustments of the solar composition. Sound speed variations at the 0.1 - 0.2 % level 
arise from reasonable changes in the input physics. This sets a threshold below which deviations 
between the Sun and the models can be explained by known observational and theoretical sources. 

It is possible to construct different relative mixtures of the light species CNONe which have 
similar opacities at the base of the convection zone, and thus similar scalar properties. Changes in 
C and N are weakly constrained because their ionization structures are similar to that of oxygen 
and they are much less abundant (by the same token, even large variations in C and N cannot 
produce a low surface oxygen). Neon, by contrast, is abundant enough that large variations could 
change the estimated oxygen, and uncertain enough that such variations cannot be dismissed 
out of hand. We find that changes in the Ne/0 level produce a strong differential sound speed 
signature atR = 0.56Rq which is explained by Ne retaining electrons (and thus higher opacity) to 
higher temperatures than O does. Although the resulting bounds are modest (0.08 dex in oxygen), 
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there would be a clear signal from a very high Ne / very low O mixture which are clearly not 
present in the data. 

Our method is sensitive to opacity, and as a result the most likely systematic error source is in 
the theoretical quantum mechanical opacity calculations despite the very good agreement between 
independent theoretical calculation. Experiments are underway to check on this possibility 
and transition to an abundance scale primarily based on experiments. However the same high 
abundances are inferred from two other independent helioseismic diagnostics - the level of 
metal ionization in the convection zone and the mean molecular weight of the solar core. In a 



comprehensive review, 



Basu & Antial ( 2008 ) found that increases in neon did not reduce the need 



for a high oxygen in the convection zone (c.f. their Figure 21). One would therefore need to 
invoke multiple errors in distinct interiors and atmospheric input physics, with the same sign and 
magnitude, to produce an acceptable high Ne - low O mixture. 

Our best interiors es timate for conipositio n is A{0) = 8.86 ± 0.04, A{Ne) = 8.15±0.17 



and A(Fe) = 7.50 ±0.05. 



Lodders et al. 



(|2009r) combined recent atmospheres measurements 
for estimates of A(C>) = 8.73 ±0.07, AiNe) = 8.05 ±0.10 and A(Fe) = 7.45 ±0.08; these two 
scales are consistent within the errors. A weighted mean of the two leads to A{0) = 8.83 ± 0.04, 
A{Ne) = 8.08 ± 0.09 and A{Fe) = 7.49 ± 0.05 which we contend is the most precise current 
estimate for abundances. In DP06 we derived relationships quantifying combinations of 
C, N, O, and Ne with the same Rcz and K,„r/- In the case of Ne, nonlinear effects can be 
induced when the change in the Ne/0 ratio is too large and a better fitting relationship is 
A{0) = -7.52 X {Ne/0)^ + 9m x (Ne/O)^ -3.10 x (Ne/Of -0.32 x (A^e/C>) + 8.98 

Our abundance scale is consistent with some photospheric abundance measurements 



( Caffauetal. 



20101; 



Pinsonneault & Delahave 



2009) but not the lower AGSS09 values. This 



is not a conflict between modem and primitive atmospheres treatments, or between one- and 
three-dimensional studies, but rather reflects differences between competing atmospheres models 
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and judgment calls on the choice of indicators, continuum levels, and the proper treatment 
of blending features. We are therefore hopeful that helioseismology may provide an absolute 
abundance standard which can be used to discriminate between competing models and which 
can be used to calibrate the appropriate composition diagnostics for the next generation of stellar 
models. 
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grant DE-FG52-09NA29580. MP is grateful to the Observatoire de Paris for their hospitality and 
support during a visit in April- June 2010. 
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Fig. 1. — Top Panel: Effect of physics changes on the predicted sound speed profile for models 
with the composition adjusted to reproduce the same Rcz and Ysurf. Detailled Uncertainties: Dot- 
Dashed (Magenta) -^ Settling ±16% , Solid (Black) -^ kr(OPvsOPAL), dotted (Red) -^ criRe), 
square (Green) -^ Age ±0.0\Gyr, Y (Blue) -^ L© ±4%, long dashed (Cyan) -^ EOS: Bottom 
panel: errors in the inverted solar sound speed profile (light gray), sound speed errors arising 
from uncertainties in the fitted Rcz and Ysurf (medium light gray), theoretical errors from the top 
panel (medium dark gray) and the total error budget (dark gray) from adding all of the above in 
quadrature. 
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Fig. 2. — Models with different Ne/O ratio but the same Rcz and Ysurf are compared. These models 
have Ne and Fe-group abundances adjusted to fit R^ and Yq . Top Panel: They are compared 
to a reference model (see text). Bottom panel: The same models are compared to the Sun. The 
dashed line corresponds to our preferred model with {0,Ne,Fe,Ne/0) = 8.86,8.15,7.50,0.197) 
which minimises the rms for the sound speed. Cases with A(0) from 8.90 (top) to 8.70 (bottom) 
in 0.05 dex increments are also shown (solid lines). 
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Fig. 3. — Published solar mixtures compared with the scalar Rcz and Ysurf constraints (top panel) 
and sound speed profile (bottom panel). All models have a settling coefficients reduced by 20% to 
account for rotational mixing. In this figure AGS05 -^ Asplund et al. 2005, AGSS09 -^ Asplund 
et al. 2010, CaffaulO -^ Caffau et al. 2010, GS98 -^ Grevesse & Sauval 1998 and Best Ne/0 is a 
modified GS98 mixture with Fe-group and Ne adjusted for seismic constraints. 



